Relativistic Effects in the Pulse Profile of the 2.5 msec X-Ray 

Pulsar SAX J1808.4-3658 
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ABSTRACT 

We analyze the properties of tlie pulsed emission from the accreting millisecond 
pulsar SAX J1808. 4-3658 in observations of its April 1998 outburst by the Rossi X-Ray 
Timing Explorer. Pulse phase spectroscopy shows that the emission evolves from a 
hard spectrum (power law with photon index 2.39 =b 0.06) to a soft spectrum (index 
3.39 ± 0.24). This softening is also observable as a phase lag in the fundamental of low- 
energy photons with respect to high-energy photons. We show that this lag is roughly 
constant over ten days of the outburst. We fit these data with a model where the pulse 
emission is from a hot spot on the rotating neutron star and the flux as a function of 
phase is calculated including the effects of general relativity. The energy-dependent lags 
are very well described by this model. The harder spectra at earlier phases (as the spot 
approaches) are the result of larger Doppler boosting factors which are important for 
this fast pulsar. 

Subject headings: accretion — black holes - stars: neutron — X-rays: stars 



1. Introduction 

Strong X-ray pulsations with a 2.49 msec 
period were discovered from SAX J1808.4- 
3658 in an April 1998 observation with the 
Rossi X-ray Timing Explorer (Wijnands &: 
van der Klis 1998). The pulsar is a member 
of an accreting binary system with an orbital 
period of 2.01 hour and a low-mass compan- 
ion (Chakrabarty & Morgan 1998). Though 
similar to other low-mass X-ray binaries in its 
timing and spectral properties (e.g. Wijnands 
k van der Klis 1999; Heindl & Smith 1998), 
SAX J1808.4-3658 is unique for its X-ray pul- 
sations. No other such binary has shown co- 



herent pulsations in its persistent flux despite 
of careful searches (Vaughan et al. 1994b, and 
references therein). 

As such, SAX J1808.4-3658 is the fastest 
rotating accreting neutron star. If the pul- 
sations are due to modulated emission from 
one hot spot on the neutron star surface, the 
2.49 msec period corresponds to an equatorial 
speed of approximately 0.1c. With these high 
speeds, SAX J1808. 4-3658 offers an excellent 
system for studying relativistic effects. 

One such effect may be the observed lag 
of low-energy photons relative to high-energy 
photons in the pulse discovered by Cui et al. 
(1998). Cui et al. (1998) suggest that the lags 
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are due to Comptonization in a relatively cool 
surrounding medium. Alternatively, the lags 
may be the result of a relativistic effect: the 
high-energy photons are preferentially emit- 
ted at earlier phases due to Doppler boosting 
along the line of sight. This possibility was 
suggested for the similar lags in the 549 Hz 
oscillations in an X-ray burst of Aql X-1 (Ford 
1999), where a simple model showed that the 
delays roughly match those expected. 

In the following, we present new mea- 
surements of the pulsed emission from SAX 
J1808.4-3658. We show that the energy- 
dependent phase lags are equivalent to a hard- 
ening pulse profile. We model this behavior 
in terms of a hot spot on the neutron star, 
including relativistic effects. 

2. Observations &; Analysis 

We have used publicly available data from 
the proportional counter array (PCA) on 
board RXTE in an 'event' mode with high 
time resolution (122|Ltsec) and high energy 
resolution (64 channels). The observations 
occurred from April 10 1998 to May 7 1998, 
when the source was in outburst. 

We generate folded lightcurves in each 
PCA channel. This is accomplished with the 
fasebin tool in FTOOLS version 4.2, which ap- 
plies all known XTE clock corrections and cor- 
rects photon arrival times to the solar system 
barycenter using the JPL DE-200 ephemeris, 
yielding a timing accuracy of several ^usec 
(much less than the phase binning used here) . 
As a check, we have applied this method to 
Crab pulsar data and the results are identical 
to Pravdo et al. (1997). To produce pulse pro- 
files in the neutron star rest frame, we use the 
SAX J1808.4-3658 orbital ephemeris found by 
Chakrabarty & Morgan (1998). An example 
folded lightcurve is shown in Figure 1 (top) 
for the observation of April 18 1998 14:05:40 
to April 19 1998 00:51:44 UTC. 



To study the energy spectra at each phase 
bin, we take the rates at pulse minimum and 
subtract it from the rest of the data at other 
phases. This effectively accomplishes back- 
ground subtraction and eliminates the un- 
pulsed emission which we do not wish to con- 
sider. Note, however, that the pulsed emis- 
sion may have some contribution even at the 
pulse minimum and this subtraction scheme 
represents only a best approximation to the 
true pulse emission. We generate detector re- 
sponse matrices appropriate to the observa- 
tion date and data mode using pcarsp v2.38, 
and use XSPEC v. 10.0 to fit model spectra. 

We fit the spectrum here with a simple 
powerlaw function. Though the function itself 
is not meant to be a physical description, the 
powerlaw index provides a good measure of 
the spectral hardness. Fits in each phase bin 
have reduced of 0.7 to 1.6. Including an in- 
terstellar absorption does not substantially af- 
fect these results. The powerlaw index clearly 
increases through the pulse phase (Figure 1, 
bottom), i.e. spectrum evolves from hard to 
soft. 

We also fit the profiles of the folded 
lightcurves in each channel using Fourier func- 
tions at the fundamental frequency and its 
harmonics. Prom these fits we determine the 
phase lag in each channel relative to the fits 
in some baseline channel range. Results for 
the 18 April 1998 observation are shown in 
Figure 2 as solid symbols. Note that nega- 
tive numbers indicate that high-energy pho- 
tons precede low-energy photons. We are also 
able to measure lags in the first harmonic, 
and find that they are opposite in sign to the 
fundamental, i.e. low-energy photons precede 
high-energy photons in the first harmonic. No 
lags are measurable in the other harmonics. 

Another way to measure energy-dependent 
phase lags is by Fourier cross-correlation anal- 
ysis. This is the method used for SAX 
J1808.4-3658 by Cui et al. (1998) and for 
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other timing signals as well, e.g. kilohertz 
QPOs (e.g. Kaaret et al. 1999). For a de- 
scription of cross correlation analysis see e.g. 
Vaughan et al. (1994a). Prom the PCA 
event mode data we calculate Fourier spec- 
tra in various channel ranges with Nyquist 
frequency of 2048 Hz and resolution of 0.25 
Hz. We then calculate cross spectra defined 
as C{j) = X^(j)X2{j), where X are the com- 
plex Fourier coefficients for two energy bands 
at the pulsar frequency vj. The phase lag be- 
tween the two energy bands is given by the 
argument of C. We measure all phase delays 
relative to the unbinned channels range 5 to 
8, i.e. 1.83 to 3.27 keV for 5 detector units in 
PCA gain epoch 3 (April 15 1996 to March 22 
1999). The results for the 18 April 1998 ob- 
servation, are shown in Figure 2 by the open 
symbols. The phase lags are consistent with 
those calculated from the lightcurve fitting. 
Prom the cross-correlation spectra we are not 
able to measure lags in the much weaker har- 
monics. 

We have calculated phase lag spectra also 
for other RXTE observations during the April 
1998 outburst. These spectra are similar to 
that in Figure 2. To quantify the trends, we 
compute an average phase delay, (pavg, over 
all energies for each observation. We also fit a 
broken powerlaw function to each phase delay 
spectrum: (p = E~'^ below a break energy, 
and (j) = 0maa; above this energy. Figure 3 
shows the quantities , and ^max versus 
the time of each observation. 

There is a clear connection between the 
results of the two analyses presented here. 
The phase resolved spectroscopy shows that 
the spectrum softens and correspondingly the 
peak of the pulse profile appears slightly ear- 
lier in phase for higher energies (Figure 1). 
The method of measuring phase delays shows 
the same behavior: higher-energy photons 
preferentially lead lower-energy photons in the 
fundamental and the magnitude of this phase 



delay increases with energy (Figure 2). In the 
following we discuss a model that can account 
for the phase delays/spectral softening mea- 
sured here. 

3. Model 

Wc calculate the expected luminosity as 
a function of phase in a manner similar to 
Pechenick et al. (1983); Strohmayer (1992) 
but including Doppler effects (Chen & Sha- 
ham 1989) and time of flight delays (Ftaclas 
et al. 1986). This treatment is based on a 
Schwarzschild metric, where the photon tra- 
jectories are completely determined by the 
compactness, R/M. The predicted luminosity 
as a function of phase from our code matches 
the results of Pechenick et al. (1983) and Chen 
k, Shaham (1989) for the various choices of pa- 
rameters. Braje et al. (2000) recently devel- 
oped a model for pulse profiles using a slightly 
different approach. 

The parameters in the model are the speed 
at the equator of the neutron star, v, the 
mass of the neutron star, M, the compactness, 
R/M^ the angular size of the cap, a, and the 
viewing angles, /3 (the angle between the rota- 
tion axis and the cap center) and 7 (rotation 
axis to line of sight). Another ingredient is 
the emission from the spot, which we take as 
isotropic and isothermal. The spectrum of en- 
ergy emitted from the spot is another impor- 
tant input. Note that if the emitted spectrum 
is power law like, the observed spectrum will 
not evolve with phase, since Doppler trans- 
formation preserves the shape (see Chen & 
Shaham 1989). The intrinsic spectrum must 
therefore have some shape that transforms to 
match the observed hardening and phase lags; 
we use a blackbody spectrum with tempera- 
ture UTq. 

A fit from the model is shown in Figure 2. 
This fit uses the following model parameters: 
R/M = 5, M = I.8M0 , kTo = 0.6 keV, 
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V = 0.1, /3 = 7 = 10°, a = 10°. To derive 
count rates, we use table models in XSPEC 
and the appropriate response files as discussed 
above. The fit for this single set of parame- 
ters is good; we find = 5.3 from the Fourier 
cross-correlation data or a reduced of 0.8 
for all the parameters fixed. 

A full exploration of the parameter space of 
the model is beyond the scope of this letter. 
We note the following trends, however. The 
magnitude of the phase lags depends sensi- 
tively on V and UTq. Larger delays result from 
higher speeds, because the pulses become in- 
creasingly asymmetric due to Doppler boost- 
ing. This asymmetry is also energy depen- 
dent, so there is a dependence on /cTq as well 
(see Chen Sz Shaham 1989). The lags also 
depend on /3 and 7, especially at higher ener- 
gies where there is the turn-over noted above. 
The phase lags are less sensitive to R/M and 
a. We have tested the assumption of isotropic 
emission and found that the phase lags de- 
pend only weakly on the beaming. 

4. Discussion 

The pulsed emission from SAX J 1808. 4- 
3658 evolves through its phase from a rela- 
tively hard to a soft spectrum, as shown by 
our phase resolved spectroscopy (Figure 1). 
This evolution can also be thought of, and 
measured as, an energy-dependent phase lag 
in the fundamental (Figure 2), i.e. higher- 
energy photons emerging earlier in phase than 
lower-energy photons. 

We have applied a model to the data which 
consists of a hot spot on the rotating neu- 
tron star under a general relativistic treat- 
ment. The dominant effect accounting for 
energy-dependent delays is Doppler boosting, 
the larger boosting factors at earlier phases 
giving harder spectra. The model fits the data 
quite well. The model also provides a stable 
mechanism for generating the phase delays. 



which meshes with the fact that the charac- 
teristic delays remain stable in time to within 
25% (Figure 3) even as there is a factor of two 
decrease in the X-ray flux, a possible tracer of 
accretion rate. As noted in Ford (1999), in 
addition to explaining the energy-dependent 
phase lags in SAX J1808.4-3658, this mech- 
anism may also account for the lags in the 
burst oscillations of Aql X-1 (Ford 1999) and 
kilohertz QPOs (Vaughan et al. 1998; Kaaret 
et al. 1999). Phase resolved spectroscopy has 
not yet been possible in these signals. 

The model offers a new means of measur- 
ing the neutron star mass and radius, a no- 
toriously difficult problem in accreting binary 
systems. The radius is directly related to v, 
the equatorial speed {v = Jlgpin^)- The fits 
also depend on the model parameters M and 
R/M, though the lag spectra are less sensitive 
to these quantities. The model could be inde- 
pendently constrained by future optical spec- 
troscopy of the companion which could pro- 
vide information on M and the viewing an- 
gles, a and (5. 
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Fig. 1. — Folded lightcurve in example low (top) and high (middle) energy bands and index of the 
power law of spectral fit (bottom). The profile is repeated twice for clarity. 
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Fig. 2.— Phase delays in the fundamental of SAX J1808.4-3658 relative to the 1.83 - 3.27 keV band. 
A negative number indicates that high-energy photons precede low-energy photons. Solid symbols 
indicate measurements from fitting the folded lightcurves, open symbols are data from Fourier cross- 
correlation analysis. The line shows the model for parameters R/M = 5, M = l.SM©, kT^ = 0.6 
keV, V = d.l c, and P = 'j = 10° as described in the text. 
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Fig. 3. — Parameters of the Fourier phase delay spectra for each observation. Plotted are average 
phase delay over all energies (top), the break energy (middle) and the maximum phase delay 
(bottom). 



8 



